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Abstract 

The interaction of light quarks and instanton liquid is analyzed at hnite density of 
quark/baryon matter and in the phase of nonzero values of diquark (colour) condensate. It 
is shown that instanton liquid perturbation produced by such an interaction results in an essen¬ 
tial increase of the critical value of quark chemical potential pLc which provokes the perceptible 
increase of quark matter density around the expected onset of the colour superconductivity 
phase. 


PACS: 11.15 Kc, 12.38-Aw 


The intriguing history of studying the structure of strongly interacting matter abounded in many 
interesting but somewhat contradictory phenomenological results up to the advent of quantum chro¬ 
modynamics (QCD). Since the time of its discovery the possible existence of quark-gluon plasma 
(QGP) phase is one of the most fundamental predictions of QGD. This phase is phenomenologically 
understood as a matter state where quarks and gluons being under extreme conditions of high tem¬ 
perature (T) and(or) quark/baryon density (driven by the corresponding chemical potential /r) are 
able to propagate relatively freely over the considerable (of course, on the hadron scale) distances. 
Even the simplest quantitative estimates of QGP physical characteristics were so verysimilar that 
they allowed to launch the process of QGP experimental study in ultrarelativistic heavy ion collisions. 
And although the hrst series of experiments at SPS GERN and AGS BNL have provided us with 
encouraging results a new generation of experiments is necessary to make convincing conclusions. 
The physical programme for new experimental research certainly needs serious quantitative analysis 
of the QGD phase diagram. 

It seemed lattice QGD could be able to resolve the problem starting even on the hrst theoretical 
principles. However, up to now this approach was effective and apparently reliable in producing the 
results mainly around the temperature axis of the (/r—T)-plane and, unfortunately, relies on the held 
theory formulation in imaginary time. Nevertheless, lattice QGD has corraborated the phenomeno¬ 
logical expectations that its phase structure is quite rich [|I|. In particular, nowadays we know that 
the chiral symmetry of initial QGD Lagrangian becomes broken and the quarks are getting to be 
conhned when the temperature decreases below a critical value. Moreover, the lattice calculations 
show us that both physical phenomena result from the corresponding phase transitions and their 
critical temperatures practically coincide. The latter observation indicates that both fundamental 
phenomena could be tightly knitted and it should heuristically play an important role to illuminate 
the conhnement mechanism. 

The features of the QGD matter at high baryon densities remains poorly understood in spite 
of the recent splash of renewed activity |0]. One of the major theoretical reasons comes from se¬ 
rious technical difficulties of lattice simulations with Monte-Garlo techniques due to the complex 
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fermion determinant in QCD at finite quark densities. From the phenomenological side the recent 
experimental data are not very informative for developing the theory in this direction. It is clear 
that present and future experiments with heavy ions at higher and higher energies should produce 
the strong interacting matter rather at higher temperature and closer to the temperature axis of 
the (/i—T)-plane. In such a situation the phenomenological treatment of studying along the /r-axis 
should basically rely on the astrophysical observations of the compact stars. In the meantime, it was 
demonstrated that the diquark pairing in colour anti-triplet attractive channel induced by the instan- 
ton interaction vertex should be much more effective than an one-gluon exchange and at sufficiently 
large quark/baryon densities and small T could lead to colour superconductivity and a sophisticated 
picture of the QCD phase diagram p. 

The present paper continues our preceding one where the light quark interactions with the 
instanton liquid (IL) at nonzero values of quark/baryon chemical potentials were investigated in the 
phase of broken chiral symmetry within the improved method of calculating the generating function 
of the IL theory [0. The new treatment of the functional integral is mainly related to take into 
account the IL back-reaction upon the quark presence which was always considered negligible. In 
fact, this effect is pretty weak and should manifest itself in the leading order of expansion in the 
effective coupling. Nevertheless, it was qualitatively argued that the interaction of IL and quarks 
could increase the quark matter density around the onset of expected colour superconductivity. As 
known [^, |P this density occurred to be surprisingly small (ug ~ 0.062 fm“^) and, hence, the diquark 
phase would come to play already at the density of nuclear matter {un — 0.45 fm~^ and quark matter 
density is taken to be three times larger than the conventional nuclear one). 

Here we obtain the quantitative estimate of the corresponding critical chemical potential Hc ana¬ 
lyzing the system of the Gorkov equations for the colour superconductor which we further calculate. 
We do not introduce any simplihcation of the initial Lagrangian trying to deal with the ’exact’ four- 
quark interaction which is generated by (anti-)instantons. Such an intention is simply compulsory 
because it is dictated by the precision of problem where the effects of the IL perturbation are expected 
to be of the same order of magnitude compared to the terms which could usually be ignored. 

Let us point out that in the IL approach the generating functional has the following factorized 
form 

Z = Zg ■ . 

Here the hrst factor provides us with information on the gluon condensate whereas the fermion factor 
Z^ serve to describe the quark practice in the instanton environment P], |0|. In what follows we 
use the notations of Ref. [Q where the dimensionless variables (motivated by the form of interquark 
interaction kernel) were introduced, for example, for the chemical potential it looks like /i —> /rp/2 
and for the momenta kip/2 ki, i = 1,... ,4 where p is the average size of pseudoparticle (PP). 
The quark determinant Z^ may be transformed with the auxiliary integration over the parameter A 
to the following form (for the colour SU{3) group with the quarks of two flavours Nf = 2): 


Z,i, ~ / dX 


exp < In 


np 

Yn 


- 1 


X 


\ r dk 

^ j y ^ ^/(^) 2 {-k - ipYf{k) + ^ r ’ 

R = 2A {4^ Li V’f)(V’r L2 4) + 2A {44 Ri V’f)(4'' R2 V’: 


( 1 ) 


where 'ipj = { 1 /^, 4 )^ / = 1)2 are the quark helds with the spinors of dehnite chirality, = 

R± t//) R± = ^ n is the IL density, = (0,/i) and N is the normalizing factor. For clarity 

we take it to be equal a unity but, in principle, it could play a role of free model parameter. This 
factor is inessential for the models with the hxed value of the packing fraction parameter np‘^ but in 
the model where it is admissible for variation the weak logarithmic dependence on N appears. The 
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factors 2 in Eq. (|l|) come with making use of the dimensionless variables. The term of four-fermion 
quark interaction V may be directly expressed with the chiral components as 

L/ ^f) = J , 

it being known that for the right hand helds the substitution L —y R should be performed. The 
kernels are dehned by the functions h,, i = 1,... ,4 and the zero modes (the solutions of the 

Dirac equation with the chemical potential /i in the PP held) 

hAki,k]ii) = ^{{k-^-iki)[(2ki + i^)f7+i{k-^-ik/^)f2] + 

Ak 

+ {k + ^ + ik4)[{2k4 + ifi)fi - i{k + fi + ik4)f2]} , 


hi{k4,k;n) = 
+ 
+ 


71 ki 

4A:2 


|(2fc - ij,){k - ii- ik4)f^ + {2k /i)(fc + n + ik^jf^R 


2 {k — fi){k — /i — ik/i) — —(/i -|- ik4)[k4 {p — f2 + 
2 {kfi){kfiik^)-^{pik4)[k4{p /^j’ , 


where k = |A;| if the spatial components of 4-vector k„ are considered and 

f± ^ Ii{z^)Ko{z^) - Io{z^)Ki{z^) ^ h{z^)Ki{z^) ^ (^^kl + {k±pf 

z 2 

with the modihed Bessel functions Jj, Ki {i = 0,1). Let us introduce the scalar function h{k 4 , k; p) re¬ 
f¬ 
lated to three-dimensional component ^*(^ 4 , f; p) = h{k 4 , k; p) ,i = 1,2,3 (we omit the arguments 

of functions hi when it does not mislead). 


= S“(p;,.)e“ Ur £/ye"”S+fe;-A*) , 


'iik 


kl jra r/'t<T (Tnc+ 


S{p; p) = {p + ip) h'^{p; p) , S^{p-, —p) =h (p; —p){p + ip)^ where it is valid for the conjugated 

* I 

function (p; —p) = h^{p; p) and e is an antisymmetric matrix with ei2 = —621 = 1. Here p^ and 
other similar designations are used for the four-vectors spanned by cx^-matrices, = (±*cr, 1), (cr 
is the three-vector of the Pauli matrices), for example, p^ = p'^a^ and f/ is a matrix of rotations in 
the colour space. The similar relations are valid for the right hand components 

{i’T ^/) = / V’/ap/p/) Rapf{pf,qf;h) , 

with the kernel 

R^(p,q-,p) = ur u}re’”'T+(q-, -^) , 

* 

where T{p; p) = {p + ip)^ h~{p; p) , T+(p; —p) =K^{p', —p){p + ip)~ . The components of matrices 
(p -|- ip)'^ and hh^{p-, p) commute because the vector-function h{p) is spanned by the vector p only. 
Then the following identities are easily understood 


T{p-,p) = S+ (p; -p) , T+ (p; -p) = S{p-,p) . 


In what follows we omit the p-dependence of matrices T, T+ as the chemical potential enters the 
matrix T being always positive and the matrix T+ being negative only. Besides, two other identities 
would also be helpful 

CT2T^(p)a2 = T+(p) , CT2 T+'^(p)cT2 = T(p) , 
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where means a transposed matrix. 

Dealing with the following averages interesting to study the diquark condensates (0) 

(?)) = D 2 TT^ S{p + q) , 

= Sfg Sap vr^ - q) Gf/^{p) , 

and making use of the effective action of Eq. (2) one may obtain the Gorkov equations similar to 
ones treated in Ref. P 


gUp) 

G„+(P) 

G„-(p) 

G„-(p) 


Fpp) - E«(p) G“ ’■(-p) = 0 , 
G“(p) - S«(p) F+'> ^(p) = 1 , 

F®(p) - E''(p) G'“‘ ^(-p) = 0 , 
G“(p) - E^p) F+'- ’■(p) = 1 , 


( 2 ) 


where 


G^ 


(P) 


n -1 


= —2 (p + means the free Green function, S^(p) = A'^ T{p) e T'^{—p), 
S^(p) = e T~^'^{—p) and A^’^ denotes a gap. The form of S-matrices results from the 

kernel structure of equations if averaging over colour orientations done (remember that we consider 
colour stochastic ensemble). In order to complete Eqs. (j^) we need the following gap equations 


eA^ = 


2A f dq 
iVc(iVc - 1) J ^ 


T+(g) F^{q) T+^(-g) - T+{-q) F^ ^{q) T+^(g) 


eA^ = 


2A 


N,{N, -1) J TT 


dq I 

~4 


Tiq) F^{q) T^\-q) - T{-q) F^ ^{q) T^\q) 


Apparently, the right hand sides of these equations are proportional e because they are the differences 
a matrix and its transposed form. Similar equations are valid for the conjugated matrices 


F+l '^(p) 

G“(p) 

F+R ^( p ) 
G^^(p) 


Go(p) 

Go(p) 

Gt{p) 

Gt{p) 


_GRL S+^(p) = 0 , 

-F^(p) S+^(p) = 1 , 

-G^^^{-p) S+^(p) = 0, 
-F^(p) S+^(p) = 1 , 


( 3 ) 


with the corresponding gap equations 


eA+^ 


2A 

iVc(iVc - 1) 



T^{-q) ^{q) T{q) 


T^{q) F+«(g) T{-q) 


eA+^ 


2A 

iVc(iVc - 1) 



T+^(-g) ^(g) T+(g) - T+^(g) F+^(g) T+(-g) 




where S+'^(p) = A+^ T+^(—p) e T+(p) and S+^(p) = A+'^ T'^(—p) e T(p). In this paper we limit 
ourselves to treating the diquark condensate only. However, as known the mixed phase of non¬ 
zero values both chiral and colour condensates could exist at Pc ~ 300 MeV realizing the transitional 
regime for the onset of colour superconducting phase. In order to bring this phase to the play the 
equation system should be extended including another averages as 


(V’/af (?)) = Sfg TT^ 6(p - q) GfF^^ip) . 
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From the Eqs. (|) and (||) we find 

G^^ip) = G+(p) + G+(p)S^(p)F+«^(p), 

F+^^(p) = G^^^{-p) G^{p) . 

With the anxiliary matrices G~^^{p) = T~^{—p) T{p), G^{p) = T{p) T~^{—p) we may rewrite 

the matrices S as 

S+^(p) = e C'+^(p) , S^(p) = G^{p) e . 

The identities for T-matrix mentioned above help to show the validity of the relations 

e G^^ \-p) = G^\p) , 

e G^ ^{-p) = G^{p) . 

And for the Green fnnction we have 

G“(p) = G+(p) + G+(p) C«(p) € G“ ’■(-p) € C*‘-(p) G+(p) . 

Combining the properties of G matrices together with the identities for free Green fnnctions 

CX2 Go'^ip) a2 = G^{p) , 

we are able to get the complete equation for calculating the function G^^ in the form 

€ G“(-p) = G„-(-p) + G„-(-p) C+yp) G“(p) C*(p) G„-(-p) . 

The vector-function h{p) being spanned by the vector p helps to conclude that all the matrices 
G+l^ commute with each other. Then searching the solution for the Green function G^^ by 
iterating one finds immediately that the Green function commutes with those matrices and finally 
obtains the equation to calculate it in the form 

G“(p) = [1 + H(p)] G+(p) + H^(p) G‘-’‘ip) , 

or 

|l-i/(p)| G“(p)=G+(p), 

where H{p) = Gq{p) G^{p) Gg (—p) G+'^(p). 

The structure of matrices H{p) allows us to conclude that their sum 

H{p) + H{-p) = a{p) 

and product 

H{p) H{-p) = j3{p) , 

are proportional to the unity matrices (it is clear by definition that a(—p) = a(p), /?(—p) = Pip))- 
Introducing notation = h^i—p) we are able to present the functions Q!(p) and Pip) in the 
compact form as 

a(p) = 4 A+^ [-4 A(p) ihg) + 2 (p^ + ^ 

where 

A(p) = (p2 pp) ihg) -2i p pig^h - h^g) . 

and a scalar product is naturally defined ihg) = hiP* together with the functions ih?) and (p^) 
squared. Getting the last term of A(p) we used the designation of scalar function h mentioned above 
and for the function P we have 

m = 16 (A«)" (a«)" {p‘ + pf (h^ . 
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As both functions a and (3 are spanned by the unit matrices the solution for the Green function 
may be given in 

i-a{p)+m • 


The gap equation then looks like 


2A f dp a{p) — 2 Pip) 
Nc{Nc-l) J TT^ A+^(l - a(p) +/3(p) 


( 4 ) 


We are interested in the solution of the form at A < 0, which is dictated by 

the symmetries of four-quark interaction kernels. Let us remember that any kernel for every quark 
sort traditionally carries the imaginary unit factor ? 0, [§. The sign choice of A relies on this fact 
but, in principle, there is an alternative A^ = A’*'^ = —A^ = —A"*"^ for A > 0, if the kernels are 
redehned. An analysis shows that the denominator of Eq. is always positive and, therefore, the 
solution of this equation does exist at pretty large A. 

The quark matter state at hnite chemical potential is dehned by the saddle point of functional Eq. 
(|^ which we treat further maintaining the hrst nonvanishing contribution which is the hgure-eight 
type diagram (see, for example, P), 


I 


^ ’ J 71^ l-a{p) + f3{p) ■ 


In our consideration this contribution to the generating functional {Z^ ~ e^) 


would occur to be 


W =-npMnA +A (y) , (E) ~ / . 


( 5 ) 


For the simplest situation of constant IL density the saddle point equation reads 

- A (E) = 0 . 

Comparing it with the gap equation Eq. (Q) one could notice the peculiar feature which is very 
practical to keep numerical calculations under control. The saddle point equation leads to the 
condition of gap independent of p. 

It was demonstrated in Refs. , 0 that the quark backreaction upon IL could be perturbatively 
estimated by studying the small variations of the IL parameters 6n and 6p around their equilibrium 
values n and p. Such variations are incorporated by the IL theory if the deformable (crumpled) 
(anti - )instantons of size p, being the function of x and z, i.e. p —*• p{x,z), are treated as the 
saturating conhgurations of the functional integral. In addition, the variations of zero modes in 
the interaction vertices of quark determinant at the transformation p —> p + 5p should be taken 
into account. Then as the output we have that for the long wave length excitations (for example, 
TT-mesons) the deformation held describes colourless scalar excitations of IL with the mass gap M 

of the order of several hundreds MeV, = L where z/ = b = ^ 

kinetic coefficient being derived within the quasi-classical approach. and Nf are the numbers of 
quark colours and havours, respectively. Our estimates give for this coefficient value of a few single 
instanton actions (3 = j^ i.e. k, ^ c f3 (with the factor c ~ 1.5 — 6 according to the ansatz 
taken for the saturating conhgurations) 0. 

Then, besides the diagrams with four legs (see, the term V in (^)) the extra diagrams with the 
scalar held attached (relatively speaking, the derivative in p of the vertex function in which the 
variation of the functions describing the zero mode 

dh ■ 

hi^hi + -^5p, i = l,...,4, (6) 
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should be performed) are generated. Due to the diquark condensate presence handling the Lagrangian 
is substantially simplihed if one restricts oneself with the leading contributions coming from the 
tadpole diagrams. Indeed, the leading contributions come from the term V and the term shown in 
Fig. 1 where two vertices are linked with the propagator of scalar held. 



Figure 1: The diagram of tadpole approach (see, the text). The solid (dashed) lines correspond to 
the fermion (scalar) held. 


Analyzing the modihed Lagrangian it is easy to see that the form of the Gorkov equations would 
be the same if everywhere S is meant as an improved one S + in which the correction is 
constructed with the modihed functions Eq. (|^). Moreover, the result for the Green function is 
formally the same if the modihed functions S are meant. In particular, the form of gap equation is 
also retained because its kernel includes again the diherences of the matrices which are constructed 
from the derivatives in p of matrices T, T+ and their transposed forms. This diherence is spanned 
by the unit vector as well. However, the following substitutions 


a{p) a{p) + 5a{p) , 


where 


(3{p) (3{p) + 5f3{p) , 

da 


5a{p) = J 


dp 


sm = j^, 

are supposed to be performed and here J is the contribution of hgure-eight type diagram with the 
scalar held propagator as an external leg 


J = 


2 A 


N^{N^-1) np^K AM- 


Let us mention that ehectively a dependence on the kinetic term k disappears (remember that 
= ■^) and its precise value is not operative in the approximation developed. The derivative of 
the hgure-eight type diagram which will shortly be necessary to proceed looks like 

dl f dp \ 1 —/3(p) a{p) — 2 

dp JM [{l-a{p)+(3{p)y [1 - a{p) + f3{p)) 

And hnally the modihcation of generating functional should be accommodated bringing Eq. (§) to 
the form 

W = -np^ l^ln ^ - 


1 + A (F) . 
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Then the new (if the variation of the IL density is absorbed) equation to calculate the saddle point 
reads 

-4 

np^ — A {np^y In — A {Y) = 0 . 

and the IL density is @] 


-4 

np = 


2pe 


u y {%) Tiu + l/2) 

2pe) pe 


(7) 


where the prime available means the derivative in A, the constant ^ tt^ is a measure of 

i 1 

interaction in the stochastic ensemble of PPs with an average size as 

pA = exp The derivative in A of the hgure-eight type diagram looks like 

A' f dp f 1 + a — 6/5 + aP + —4 + 7a + 4/? — aP — 1 

OCX ^- 0(j / 




(1 — a + Pp 

and the derivative A' is dehned by the following equation 

NPN^ - 1) 2A' [dp 2/3 - 2/32 + 2ap - 


(1 — a + py 


2A2 


A^ 7 TT^ (1 — a + py 



Figure 2: The saddle point A as the function of chemical potential p at = 3 and Nf = 2. The solid 
lines correspond to the phase of broken chiral symmetry (the lower one absorbs the IL perturbation). 
The dashed lines correspond to the colour superconductivity phase (the left hand one does not include 
the tadpole contribution whereas the right hand one does. 


Fig. 2 shows the results of calculating the parameter A (which is proportional to the free energy 
within the precision of one loop approximation) as the function of chemical potential at Ac = 3, A/ = 
2. The dashed lines expose a behaviour in the phase of nonzero values of diquark condensate (the 
left hand line (lower line) corresponds to the calculations if the quark interaction with IL is ignored) 
whereas the solid lines show the behaviours in the phase of broken chiral symmetry (upper line for the 
quark interaction with IL ignored) Q. The saddle point parameter Ai of Ref. [§] and that exploited 


in the present paper are related as A^ = 


A np 


^4 


The crossing points of the curves are hxing 


2(Ac - l)Ac 

the onset of colour superconducting phase. As mentioned above, actually, there is a transitional 


















region of mixed phase where both chiral and diquark condensates develop non-zero magnitudes and 
the descent is not so steep. However, these details are inessential for further discussion. What is more 
interesting to be noticed comes from the crossing point of lower solid and dashed lines positioned on 
the plot. 



Figure 3: The quark density being dehned by the crossing point of the upper solid and left hand 
dashed lines in Fig. 2. The solid line corresponds to the stable phase and the dashed one corresponds 
to the metastable one. 

The hrst crossing point (along the /x-axis) which corresponds to the approach without the tadpole 
contributions included looks to be slightly larger in magnitude than that obtained earlier /Xc — 
300 MeV in Refs. [§], [^]. However, it does not signal about shifting the critical point to larger values 
of /X in the approach dealing with the initial Lagrangian. The IL parameters of various approaches 
are sometimes slightly different and, in principle, could be optimized htting, for example, Aqcd- 
Thus, the shift of /Xc value pointed out is well within a precision of the IL theory. Regarding the 
second crossing point (when IL is perturbed by quarks) the conclusion could be more indicative. 
Here /x becomes significantly larger. The corresponding quark matter densities as the functions of 
chemical potential for both approaches are depicted in Fig. 3 and Fig. 4. Apparently, on both plots 
the density values corresponding to the onset of colour superconductivity phase are noticeably larger 
than the density of normal nuclear matter. Perhaps, our result could be taken as a general indication 
on the considerable role of supplementary (insignihcant on the instanton background) interactions of 
light quarks. They are able to change the estimate of critical density for colour diquark condensation 
drastically. 

Now we would like to comment on the crossing point of upper dashed and lower solid lines which 
is not shown in Fig. 2. It corresponds the unrealistic values of /x and in that region the IL is hardly 
applicable. One of the limits for IL approach coming from the very large values of chemical potential 
points out that with quark matter density increasing the average interquark distances may become 
very small and the ’Coulomb’ (perturbative) held strengths would occur to be comparable with 
the (anti-)instanton ones. Then the (anti-)instanton superposition is not a proper conhguration to 
saturate the functional integral. Thus, our conclusion resulting from the approach developed gives a 
message that the quark perturbation of IL leads the corresponding curves to move apart inherent in 
the phases of non-zero values of diquark and chiral condensates in Fig. 2. The ’chiral’ curve becomes 
steeper but ’diquark’ curve is displaced to the larger values of /x increasing its critical value. It looks 
like giving more reliability to our understanding of the perturbative helds role in the IL theory could 
provide us with more accurate estimates for the phenomena considered. 

Fig. 5 demonstrates the IL density as the function of /x when the quark backreaction is incor¬ 
porated in the phase of broken chiral symmetry (solid curve) and in the phase of non-zero values of 
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Figure 4: The quark density being defined by the crossing point of the lower solid and left hand dashed 
lines in Fig. 2. The solid line corresponds to the stable phase and the dashed one corresponds to the 
metastable one. 


2.0 



1.0 


0.5 



0 400 800 '200 


jd MeV 


Figure 5: The IL density for the approach at iVc = 3 and Nf = 2. The solid line corresponds to the 
calculation in the phase of broken chiral symmetry. Both dashed lines correspond to the calculation 
in the phase of non-zero diquark condensate values. The lower dashed line gives the estimate of IL 
density without the tadpole contributions included. Regarding the upper dashed line see the text. 


diquark condensate (dashed curve). It is interesting to notice that the curves manifest the different 
character of quark influence on IL. In the phase of broken chiral symmetry the density is almost con¬ 
stant in the suitable interval of fi whereas in the colour superconducting phase quickly disappears. 


51 

J-p 


m 


It is interesting that for the former the corresponding analogue of the tadpole contribution 

Eq. d^) is strictly positive and, therefore, can lead to the increase of the IL density (gluon conden¬ 
sate) which is truely insignificant and demonstrates simply the approach sensitivity to the dynamical 
quark mass variation |Q. In the latter case the sign of tadpole contribution is changing. As seen 
from Fig. 5 it leads to the drastic change of the IL density behaviour and the gluon condensate is 
getting weaker. Strictly speaking this effect has already been discussed in [|T^ 
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